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can be represented as the theories of real (complex) F(R) gravity through the inverse (complex) 
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for the obtained solutions of the crossing of the phantom divide under a quantum correction coming 
from conformal anomaly. 
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I. INTRODUCTION 



It is observationally confirmed that the current expansion of the universe is accelerating [l[ 

aa 

. Various scenarios to explain the current accelerated expansion of the universe have 



nproposed. T 



beenproposed. The mechanism, however, is not well understood yet (for recent reviews, 
see 

There are two approaches to account for the current accelerated expansion of the universe. 
One is to introduce some unknown matter, which is called "dark energy" in the framework 
of general relativity. The other is to modify the gravitational theory, e.g., in simplest case 
to study the action described by an arbitrary function of the scalar curvature R, which is 
called gravity". Here, F(R) is an arbitrary function of the scalar curvature R (for 

reviews, see 7|, |8j). 

According to the recent various observational data including the Type la supernovae Gold 
dataset ({]], there exists the possibility that the effective equation of state (EoS) parameter, 
which is the ratio of the effective pressure of the universe to the effective energy density 
of it, evolves from larger than —1 (non-phantom phase) to less than —1 (phantom one, in 

, crosses —1 (the phantom divide) 

currently or in near future. 

A number of attempts to realize the crossing of the phantom divide have been made 
in the framework of general relativity: For instance, scalar-tensor theories with the non- 
minimal gravitational coupling between a scalar field and the scalar curvature [ij] or that 
between a scalar field and the Gauss-Bonnet term one scalar field model with non- 
linear kinetic terms JjiJ or a non-linear higher-derivative one [jjj] , phantom coupled to dark 
matter with an appropriate coupling , the thermodynamical inhomogeneous dark 



model [3], multiple kinetic k-essence [l7j], multi- field models (two scalar fields model 



energy 



18 



19 



201 ]. "quintom" consisting of phantom and canonical scalar fields 2110, and the description of 



those models through the Parameterized Post-Friedmann app roach 22|, or a classical Dirac 



field 23J or string-inspired models 24J , non-local gravity 25|, |26|] , a model in loop quantum 



cosmology 27] and a general consideration of the crossing of the phantom divide 



(for a detailed review, see {fl). In fact, however, explicit models of modified gravity realizing 
the crossing of the phantom divide have hardly been investigated, although there were 
suggestive and interesting related works ,|32|. 
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In the present paper, we study a crossing of the phantom divide in modified gravity. We 
reconstruct an explicit model of modified gravity in which a crossing of the phantom divide 
can be realized by using the reconstruction method proposed in Ref. 33|. Furthermore, we 
investigate the corresponding scalar field theory, in which there exist the Einstein-Hilbert 
action and a scalar field, obtained through a conformal transformation of the modified 
gravitational theory, and compare the evolution of the universe in the modified gravitational 



theory with that in 
Big Rip singularity 



;he corresponding scalar field theory. It is shown that the (finite-time) 



34 



35| appears in the reconstructed model of modified gravity (i.e., in 



the so-called Jordan frame), whereas that in the corresponding scalar field theory obtained 
through the conformal transformation (i.e., in the so-called Einstein frame) the singularity 
becomes the infinite-time one. Moreover, we consider the relations between the scalar field 
theories with realizing a crossing of the phantom divide and the corresponding theories of 
modified gravity by using the inverse conformal transformation of scalar field theories. It 
is demonstrated that the scalar field theories describing the non-phantom phase (phantom 
one with the Big Rip singularity) can be represented as the theories of real (complex) F(R) 
gravity through the inverse (complex) conformal transformation. On the other hand, a 
very realistic model of modified gravity that evades solar-system tests and realizes a viable 



cosmic expansion in the past has recently been proposed in Ref. 36| (for some related models, 



sec 



33, 



381]) . In this model, our universe is asymptotically de Sitter space. Therefore, we also 



reconstruct a model of modified gravity in which the transition from the de Sitter universe 
to the phantom phase can occur in such a viable theory. In addition, we explore the stability 
for the obtained solutions of the crossing of the phantom divide under a quantum correction, 
in particular conformal anomaly. 

Our goal in this paper is to show that in principle the crossing of the phantom divide 
can be realized in the framework of modified gravity without introducing any extra scalar 
components with the wrong kinetic sign (phantom). We reconstruct such an explicit model 
of modified gravity. By presenting it, it can be illustrated that the crossing of the phantom 
divide can occur in modified gravity as the scalar field theories in the framework of general 
relativity. The demonstration in this work can be interpreted as a meaningful step to 
construct a more realistic model of modified gravity, which could correctly describe the 
expansion history of the universe. 

This paper is organized as follows. In Sec. II we explain the reconstruction method of 
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modified gravity proposed in Ref. [33|. Using this method, we reconstruct an explicit model 
of modified gravity in which a crossing of the phantom divide can be realized. In particular, 
we show that the Big Rip singularity appears in this modified gravitational theory. In Sec. Ill 
we consider the corresponding scalar field theory, which is obtained by making the conformal 
transformation of the modified gravitational theory with realizing a crossing of the phantom 
divide. We demonstrate that the Big Rip singularity does not appear in the corresponding 
scalar field theory. In Sec. IV we investigate the relations between scalar field theories and the 
corresponding modified gravitational ones. In Sec. V we study the viable model of modified 
gravity in which the transition from the de Sitter universe to the phantom phase can occur. 
In Sec. VI, we examine the stability for the obtained solutions of the phantom crossing 
under a quantum correction coming from conformal anomaly. Finally, some summaries and 
outlooks are given in Sec. VII. Detailed derivations and explanations about each section 
are shown in Appendixes A-F. We use units in which k-Q = c = H = 1 and denote the 
gravitational constant 8ttG by k 2 , so that k 2 = 8n/Mp 2 , where M P i = G~ 1 ^ 2 = 1.2xl0 19 GeV 
is the Planck mass. 

II. RECONSTRUCTION OF MODIFIED GRAVITY 

We investigate modified gravity with realizing a crossing of the phantom divide by using 
the reconstruction method. (The equivalence between F(R) gravity and the scalar tensor 



theory was explicitly shown in Ref. 



A. Reconstruction method 



391 ] . The limited case was given in Ref. 401 ] . 



First, we briefly review the reconstruction method of modified gravity proposed in 
Ref. [33]. 



The action of F(R) gravity with general matter is given by 

'F(R) 



S = Id rcy— g 



T -^matter 

2k z 



(2-1) 



where g is the determinant of the metric tensor and £ ma tter is the matter Lagrangian. 
The action ( 12.11) can be rewritten to the following form by using proper functions P{4>) 



4 



and Q{4>) of a scalar field 0: 



S = J d 4 X^^[P((f>)R + £ ma tter} • 



(2.2) 



The scalar field may be regarded as an auxiliary scalar field because <fi has no kinetic term. 
It follows from the action (12. ip that the equation of motion of is given by 

= » + ^, (2.3) 

which may be solved with respect to <fi as <fi = 4>{R)- Substituting <fi = <fi(R) into the action 
( 12.21) . we find that the expression of F(R) in the action of F(R) gravity in Eq. ( 12.1 ft is given 
by 

F(R) = P{<j>(R))R + Q(<f>(R)). (2.4) 
From the action (12.21) . we find that the field equation of modified gravity is given by 



-g, v [P(<P)R + Q(0)] - R^P{4>) - g^UP{<f>) + V,V,P(0) + k 2 T^^ = , (2.5) 



where V M is the covariant derivative operator associated with g^ vi □ = g^V^Vj, is the 
covariant d'Alembertian for a scalar field, and 7 i j™ attcr ) j s fa e contribution to the matter 
energy-momentum tensor. 

We assume the flat Friedmann- Robert son- Walker (FRW) space-time with the metric, 

ds 2 = -dt 2 + a 2 (t)dx 2 , (2.6) 

where a(t) is the scale factor. 

In this background, the (/i, v) = (0, 0) component and the trace part of the (/i, v) = 
component of Eq. ( 12.51) . where i and j run from 1 to 3, read 

- 6H 2 P(4>(t)) - Q(4>(t j) - 6H dP ^^ + 1k 2 p = , (2.7) 

and 

2 *p<m + 4H <msi + + <&) pw)) + <nm + = o , (2.8) 

respectively, where H = a/a is the Hubble parameter and a dot denotes a time derivative, 
= <9/<9t. Here, p and p are the sum of the energy density and pressure of matters with a 
constant EoS parameter Wi, respectively, where % denotes some component of the matters. 



Eliminating Q(4>) from Eqs. (12. 7p and (12.81) . we obtain 



d 2 P(Ht)) _ H dPW)) 



2HP(<j)(t)) + K 2 (p + p ) = 0. 



(2.9) 



dt 2 " dt 

We note that the scalar field <fi may be taken as = t because <fi can be redefined properly. 
We now consider that a(t) is described as 



a(t) = aexp (g(t)) 



(2.10) 



where a is a constant and g(t) is a proper function. In this case, Eq. (I2.9p is reduced to 
rf 2 P(0) dg{</>)dP{</>) , d 2 ~g&) 



+ 2- 



-P(<f>) 



+ k 2 (1 + ™i) Pid^ 3{1+W ^ exp [-3 (1 + Wi ) g(4>)} = 



(2.H) 



where is a constant and we have used if = dg(<p)/ (d<j)). Moreover, it follows from Eq. (I2.7P 
that Q((p) is given by 



-6 



dg{4>) 



P(<f>) - 6 



dm dp^) 



+ 2k 2 Pia' Hl+m) exp [-3 (1 + Wi ) g{<j>)\ . 



(2.12) 



Hence, if we obtain the solution of Eq. (12. lip with respect to P(4>), then we can find Q{4>). 
In Appendix A, some points on the reconstruction method are noted. 

We mention that the convenient reconstruction for scalar field theories could be given 



in Refs. 



41 



42J (for a recent review, see 43|]). Furthermore, the reconstruction in the 



scalar-Einstein-Gauss-Bonnet theories was considered in Ref. 44] 



B. Explicit model with realizing a crossing of the phantom divide 



Next, using the reconstruction method explained in the preceding subsection, we recon- 
struct an explicit model in which a crossing of the phantom divide can be realized. 
A solution of Eq. (12. lip without matter can be given by 



P(0) 



£(0) = -lOm 



tr 



C 



tr 



7+1 



p{4>) = p+<p p+ +p-<p p - 



± 



1 ± y/1 + 1007(7 + 1) 



(2.13) 
(2.14) 

(2.15) 
(2.16) 



where 7 and C are positive constants, to is the present time, and p± are arbitrary constants. 
The derivation of this solution is shown in Appendix B. 

From Eq. (I2.14p . we find that g{<fi) diverges at finite <fi when 



(2.17) 



34 



35[. (Other kinds of 



45|.) One only needs to consider 



which tells that there could be the Big Rip singularity a t t = t s 
finite-time future singularities have been studied in Ref. 
the period < t < t s because g((f>) should be real number. Eq. ( 12. 14ft also gives the following 
Hubble rate H(t): 



W d(f> \t 



7 



to 



+ (7 + 1 c 



to 



c 



7+1 



(2.18) 



where it is taken = t. 

In the FRW background (12.61) . even for modified gravity described by the action (12.11) . 
the effective energy-density and pressure of the universe are given by p c g = 3H 2 /k 2 and 
Pes = — (2H + 3H 2 j /k 2 , respectively. The effective EoS parameter w c s = p c s/p c s is 
defined as 

2H 



W e ft = -1 



3H 2 



For the case of H(t) in Eq. (12. 18ft . from Eq. (12.191) we find that w e s is expressed as 



w eff = -1 + U(t) 



(2.19) 



(2.20) 



where 



-7 + 47(7 + 1) (i) 27+1 + (7 + l)(^) 



U(t) 



2H 



2(2 7 +l) 



3H 2 




/ 1 \ 27+1" 


2 




15 


7 + (7+l) (£) 





(2.21) 



Moreover, the scalar curvature is given by R = 6 ( H + 2H 2 ) . For the case of Eq. (I2.18p . R 
is described as 



60 



R 



7 (20 7 - l) + 44 7 (7 + 1) (^) 27+1 + (7 + 1) (20 7 + 21) (£) 



2(27+1) 



t 2 



27+1 



. (2.22) 
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In deriving Eqs. (12~2T1) and (T2T22]) . we have used Eq. (l2~T7j) . 
When i — >• 0, i.e., t <C i s , behaves as 

IO7 



Hit) 



t 



(2.23) 



In this limit, it follows from Eq. ( 12.191) that the effective EoS parameter is given by 

1 

W e ff = -1 + 



157 



(2.24) 



This behavior is identical with that in the Einstein gravity with matter whose EoS parameter 
is greater than —1. 

On the other hand, when t t s , we find 

10 



H(t) 



(2.25) 



t a -t 

In this case, the scale factor is given by a(t) ~ a (t s — t)~~ 10 . When t — > t s , therefore, a — > 00, 
namely, the Big Rip singularity appears. In this limit, the effective EoS parameter is given 
by 

1 16 

w cS = -l = . (2.26) 

15 15 V ' 

This behavior is identical with the case in which there is a phantom matter with its EoS 
parameter being smaller than — 1. Thus, we have obtained an explicit model showing a 
crossing of the phantom divide. 

It follows from Eq. (12 . 1 9[) that the effective EoS parameter w c s becomes —1 when H = 0. 
Solving w c ff = — 1 with respect to t by using Eq. (I2.20p . namely, U(t) = 0, we find that the 
effective EoS parameter crosses the phantom divide at t — t c given by 

1/(27+1) 



t c = t s [ -2 7 + A /4 7 2 + 



1 



7 



(2.27) 



From Eq. (12.211) . we see that when t < t c , U(t) > because 7 > 0. Moreover, the time 
derivative of U(t) is given by 



dU(t) 
dt 



2 7 ( 7 + l)(2 7 + iy 



15 



7 + (7 + l) (i) 



27+1 



27 



l - [ T, 



27+I 



(2.28) 



Eq. (I2.28P tells that the relation dll (t) / (dt) < is always satisfied because we only consider 
the period < t < t s as mentioned above. This means that U(t) decreases monotonously. 



S 



Thus, the value of U(t) evolves from positive to negative. From Eq. (12.201) . we see that 
the value of w c s crosses —1. Once the universe enters the phantom phase, it stays in this 
phase, namely, the value of w e ff remains less than —1, and finally the Big Rip singularity 
appears because U(t) decreases monotonically. Note that other types of the finite-time 
future singularities in modified gravity are possible as demonstrated in Ref. 46]. 
It follows from Eqs. (12331 . (l234j) . (12351) and (12371) that P(t) is given by 



P(t) 



Using Eqs. ([232]) and (l2T29|) . one gets 



ts) 



(2.29) 



Q(t) 



-6H 



\ 27+1 



E(> + 7 

j=± ^ 



(2.30) 



If we can solve Eq. ( 12.221) with respect to t as t = t(R), in principle we can obtain the 
form of F(R) by using this solution and Eqs. (12.41) . (12.291) and (I2.30p . In fact, however, 
for the general case it is difficult to solve Eq. (I2.22p as t = t(R). Hence, as an solvable 
example, we show the behavior of t 2 F(R) as a function of R = t 2 s R in Fig. 1 for 7 = 1/2, 
p + = — p_ = 0, j3 + = (l + 2^/l§) 1 2 and t s = 2t . The quantities in Fig. 1 are shown 

in dimensionless quantities. The horizontal and vertical axes show R and t 2 F, respectively. 
(Here, R = t 2 s R = 4R/Rq, where Rq is the current curvature. In deriving this relation, we 
have used t s = 2to, to ~ Hq , where Hq is the present Hubble parameter.) From Fig. 1, we 
see that the value of F(R) increases as that of R becomes larger. 

To examine the analytic form of F(R) for the general case, we investigate the behavior 
of F(R) in the limits t -> and t -> t s . When t -> 0, from Eq. we find 



6O7 (2O7 - 1) 
R 



(2.31) 



In this limit, it follows from Eqs. (|2~23|) . (|2T29|) . (l23jQ|) and (l2~3ll that the form of 

F(R) is given by 



F(R) 



W6O7 (20 7 - 1)R~ 1/2 



}- s ^607 (2O7 - 1)R^/ 2 
5 7 - 1 - f3j 



27+1 



R 



3=± 



207 



Pj- [60 7 (20 7 - 1)]^ /2 i?"^ /2 



(2.32) 




FIG. 1: Behavior of t 2 s F(R) as a function of R for 7 = 1/2, p+ = -l/tf + , p_ = 0, (3+ 
(l + 2-V/19) /2 and t s = 2t . 



On the other hand, when t — > t s , from Eq. (12.251) we obtain 



t~t,-3 



140 



'2.33) 



In this limit, it follows from Eqs. Q, d2T25|) . (TO?]) . (1230]) and fT2~33l that the form of 
F(R) is given by 



F(R) 



{±[t s -3VU0R^]y 


1 - 


[l 3^40^-1/2" 


27+1 



V 



'20 



J 



R^pj [t s - 3VU0R~ 1/2 3 



j'=± 



gt a + 0%- 15) 7T 1 / 2 



1 



t s - 3VUOR- 1 / 2 



[2.34) 



The above modified gravity may be considered as some approximated form of more realistic 
viable theory. For large R, namely, t 2 s R ^> 1, the expression of F(R) in ( 12.341) can be 
approximately written as 

5 



F(R) « - 



1 



7 L3V140(27 + 1) V*o 



t. 



VJ=± 



'2.35) 
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III. CORRESPONDING SCALAR FIELD THEORY 



In this section, motivated by the discussion in Ref. [47j, we consider the corresponding 
scalar field theory to modified gravity with realizing a crossing of the phantom divide, which 
is obtained by making the conformal transformation of the modified gravitational theory. 

By introducing two scalar fields ( and £, we can rewrite the action (12. ip to the following 
form 

S = J d 4 X^ {^2 K ( R ~ + F(0\ + Matter} • (3.1) 

This is the action in the Jordan-frame, in which there exists a non-minimal coupling between 
£ and the scalar curvature R. The form in Eq. (13. ip is reduced to the original one in Eq. (12. ip 
by using the equation ( = R, which is the equation of motion of one auxiliary field £. 
We make the following conformal transformation: 

SV 9v» = eCT SV > ( 3 - 2 ) 

where 

e° = F'(C) • (3.3) 

Here, a is a scalar field and a hat denotes quantities in the Einstein frame, in which the 
non-minimal coupling between £ and R in the action (13.11) disappears. 

By defining ip as <p = a/3/2o"/k, we obtain the following canonical scalar field theory: 



5" S T = I d x 



R I 

2^2 ~ 2 



- - -r v d^d v y - V(<f) + e-V 2 /3^£ mattcr 



(3.4) 



The detailed derivation of the action (13. 4p is given in Appendix C. 
We now investigate the case in which F(R) is given by 

F(R) = ciM 2 (JL^ , (3.5) 

where c\ is a dimensionless constant and M denotes a mass scale. The form of F[R) in 
Eq. (12.351) corresponds to the one in Eq. (13. 5p with n = —7/2. It may seem that such a 
model may have problems in the description of the past universe evolution. However, there 
is the trick to make its past evolution consistent with observations described in Ref. 33] . 



It uses the introduction of compensating dark energy dominated at intermediate universe 

11 



which disappears effectively at current universe. Our primary purpose in this work is current 
universe admitting the phantom divide crossing in modified gravity, so we will not discuss 
the past evolution of the model under discussion. In this case, the scale factor a(t) and the 



scalar curvature R are given by 



47) 



and 



R 



a(t)=a(t s -tt +l){2n+1)/i ' n+2) 



6n{n + l)(2n + l)(4n + 5) 1 



(3.6) 



(n + 2) 2 

respectively. From Eqs. (13.21) and (13. 3p . we find 



(t, - t) 



2 • 



(3.7) 



dt = ±e C7/2 dt 
e -/ 2 = v /3^ 



(n + 2) ; 



(n+l)/2 



M n+i (t s - t) 



n+l 



(3.8) 
(3.9) 



_6n(n + l)(2n + l)(4n + 5)_ 

where we have used Eq. (13.71) . It follow from Eq. (13.91) that the relation between the cosmic 
time in the Einstein frame t and that in the Jordan frame is given by 



t = T 



n + 2 



(n + 2)' 



(n+l)/2 



M n+l {t a - 1) 



n+2 



(3.10) 



_6n(n + l)(2n+ l)(4ra + 5)_ 

If n < —2, the limit of t — > t s corresponds to that of t — > =Foo. For the case of Eq. (12.351) . 
n = —7/2. From Eqs. (13. 2h and (13.31) . we also find that the metric in the Einstein frame is 
expressed as 



ds 2 = e a ds 2 = -dP + a (i) dx 2 , 
where a (£) is the scale factor in the scalar field theory given by 



(3.11) 



(i) = at 



3[(n+l)/(n+2)] 2 



(3.12) 



a = a =F 



1 



n + 2 



-3[(n+l)/(n+2)] 2 



x < V—nci 



{n + 2f 



_6n(n + l)(2n + l)(4n + 5) 



(n+l)/2 



(2n 2 +2n-l)/(n+2) 2 



M 



n+l 



. (3.13) 



For n = —7/2, because when t — > t S} t — > =Foo, it follows from Eq. (I3.12p that the scale 
factor in the scalar field theory a (t) diverges at infinite time. 
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Consequently, the 'finite-time' Big Rip singularity in F(R) gravity, i.e., in the Jordan 
frame, becomes the 'infinite-time' one in the corresponding scalar field theory obtained 
through the conformal transformation of the theory of F(R) gravity, namely, in the Einstein 
frame. This shows the physical difference of late-time cosmological evolutions between the 
theory of F(R) gravity and the corresponding scalar field theory, which are mathematically 
equivalent theories. 



IV. RELATIONS BETWEEN SCALAR FIELD THEORIES AND THE CORRE- 
SPONDING THEORIES OF F(R) GRAVITY 



In this section, following the considerations in Refs. 47|, |49] , we investigate the relations 
between scalar field theories and the corresponding transformations to F(R) gravity. 
The action of scalar field theories in the Einstein frame is given by 



d x\/—g 



7^ T ln X d„X ~ W(X) 



(4.1) 



Here, in the non-phantom phase the sign of the kinetic term is — , while in the phantom one 
that is +. 

To study the corresponding theories of F(R) gravity, we make the inverse conformal 
transformation of the action of the scalar field theories (j4.ip . In the non-phantom phase, we 



use the inverse conformal transformation 



49j in order to vanish the kinetic term of x 



(4.2) 



As a consequence, the action in the Jordan frame for the non-phantom phase is given by 

-F NP (R) 



5*, 



NP 



d Xy/—g- 



2^ 



F NP (R) = e ± V^x( R ) R _2K 2 e ±2 V^x(R) ] y(x(R)) . 



(4.3) 
(4.4) 



The scalar field x is j us t an auxiliary field and can be expressed in terms of the scalar 
curvature as x — x(F) by solving the equation of motion of x' 



(4.5) 
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Similarly, in the phantom phase we use the complex conformal transformation 47| in 
order to vanish the kinetic term of x 



9r - 9^ = ^ lliRX 9^- (4-6) 
As a result, the action in the Jordan frame for the phantom phase is given by 



S P = I d\^rg-^± , (4.7) 



-FAR) 

2k 2 

The equation of motion of x is given by 

R = e^^x ^Ak 2 W(x) T iVGk^^J . (4.9) 

This equation can be solved with respect to x as X — x{R). 

In general, scalar field theories describing the non-phantom (phantom) phase can be 



represented as the theories o 



conformal transformation 



47 



real (complex) F{R) gravity through the inverse (complex) 



49j. We note that the consideration of this section can be 
applied to not only the model in Eq. (12.181) but also to any other scalar field theories 
with/without the crossing of the phantom divide: e.g., the case in which the Hubble rate 
is given by Eq. (IB4I) . In Appendix D, we examine the more detailed relation between the 
scalar field theories with realizing a crossing of the phantom divide and the corresponding 
theories of F(R) gravity. 



V. MODEL OF F{R) GRAVITY WITH THE TRANSITION FROM THE DE SIT- 
TER UNIVERSE TO THE PHANTOM PHASE 

In this section, we reconstruct a model of F(R) gravity in which the transition from the 
de Sitter universe to the phantom phase can occur by using the method explained in Sec. II 
A. 



A. Reconstruction of the viable F(R) gravity 



The interesting viable model is proposed in Ref. 36(. It is known that the above model 
is a very realistic modified gravitational theory that evade solar-system tests, which was 

11 



mentioned also in Ref. 50] . As shown in the above reference 36j, this model could reproduce 
the viable cosmic expansion, correctly describing the phases before dark energy epoch. Hence 
our universe is asymptotically de Sitter space. The form of the model in Ref. 36J and its 
generalization is presented in Appendix E. 

As an example realizing the transition from the de Sitter universe to the phantom phase, 
we can consider the following form of the Hubble rate: 



H = g 



9i 



U -t 



(5.1) 



where go, gi and t s are positive constants. When t — > — oo, H goes to a constant H — > go. 
Hence the universe is asymptotically de Sitter space. On the other hand, when t — > t s , the 
second term on the right-hand side of Eq. (15.11) dominates and H behaves as H ~ gi / (t s — t). 
It follows from H ~ g\j (t s — t) 2 > and Eq. (12.241) that w c r < — 1, namely, the universe 
enters the phantom phase. Then there appears the Big Rip singularity at t — t s . For the 
case of Eq. (15.11) . R is given by 



R = 6 



4#o#i , 9x (2#i + 1) 



(5.2) 



(t s -ty 

We consider the case in which the contribution from matter could be neglected. We take 



into account it later. Eq. (15. ip shows 

dg(4>) 



go + 



gi 



(5.3) 



where we have taken <fi = t. Substituting Eq. (I5.3P into Eq. (12.111) . we obtain 



d 2 P{(f)) 



go 



gi 



dP{<t>) 



2<7i 



lP{<t>) ■ 



(5.4) 



The solution is given by 



P{z) = C + z a F K (a, f, z) + C-Z^Fk (a - 7 + 1, 2 - 7; z) 



(5.5) 



where 



z = go 



a 



l- 9l ± y/g\ - 10 5l + 1 



7 = 1 ± 



y/gf - Wx + 1 



F K (a, 7; z) = — 



a(a + 1) • • • (a + n — 1) z n 



n=0 



7(7+ 1) • ■ ■ (7 + n — 1) n! 



(5.6) 



15 



Here, Fk is the Kummer functions (confluent hypergeometric function), and C+ and C_ are 
dimensionless constants. Using Eqs. (12.121) and (]5.5p . we obtain 



Q(z) = -6g 2 (l - |) I C + z«- 1 



a 



G-z 



(a- g 1 + z) F K (a, 7; z) + -Z-Fr: (a + 1, 7 + 1; z) 

7 

;i - 7 - g-i + z) F K (a - 7 + 1, 2 - 7; z) 



OL — 'Y + 1 

+ r , ^(a-7 + 2,3-7;z 
2-7 



(5.7) 



In Eqs. (15.51) and (15. 7p . because we take 4> = t, z = g (t — t s ). 

It follows from Eq. (15.21) that when t — ► — 00, namely, in de Sitter phase, R becomes 
constant as R ~ 12<7q. On the other hand, when i — > t s , i? ~ 6(71 (2^ + 1) / (t s — t) 2 . Using 
this relation and z — go (t — t s ), we find 



-go 



6<7i (2 9l + I] 



R 



(5.8) 



In the limit of t — > t s , \z\ <C 1 because i? diverges. Expanding the Kummer functions in 
Eqs. (15 .5j) and (15.71) and taking the first leading order in z, from Eqs. ( 12. 4ft and ( 15.81) we find 
that the form of F(R) in this limit is approximately expressed as 
R 



F(R) 



(7 - 1) {2 9l + 1) 



C+{a-l)(a + g 1 + l) -g ov / 6g 1 (2 9l + 1 



R -a/2 



C_( a -j)(2-j + 9l ) -g y/6 gi (2 9l + l) R^~^ 2 \ . (5.9) 



1-7 



Next, we study the case in which there exists the matter. We consider the cold dark 
matter with w = 0. We numerically solve Eq. (I2.1ip with the cold dark matter. To execute 
this, for simplicity, we set t s g = gi in Eq. (15. ip . We show the behavior of F(R)/ (2k 2 ) in 
Fig. 2. From Fig. 2, we see that F(R) increases in terms of R. The detailed explanation of 
the numerical calculations is given in Appendix F. 



B. Addition of a non-local term to the viable modified gravity models 

It seems very difficult in the framework of F(R) gravity to construct a model generating 
the transition from de Sitter space to the phantom phase in the viable model of the previous 
sub-section. One may add any term, which is a function of R, to such models. The term 
should be small in the present universe and may be small even in the past universe, where 
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FIG. 2: Behavior of F(R)/ (2k 2 ) as a function of R. Legend is the same as Fig. 3. 

the curvature could be large. Hence the term should dominate only at small curvature, that 
is, smaller one than the present curvature. In the asymptotically-de Sitter model above, 
however, the present universe is asymptotic de Sitter space, where the curvature is (almost) 
constant. Thus, the curvature cannot become smaller than the order of the curvature in the 
present universe and the added term never dominates. 

Let us add a non-local term to any model from Ref. {stJ by using a proper function K: 

5f — K (—\JT l R) R m , (5.10) 

where m is a positive constant. In the de Sitter universe where the curvature and the Hubble 
rate are constant R = Ro, H = Ho, we find 

D-'R = -^-t + c 4 e~ 3Hot + c 5 , (5.11) 

where c 4 and c 5 are constants of the integration. For large t, we find — D~ 1 R ~ [R / (3H )} t. 
If we choose K to be a slowly increasing function, 5f dominates in the future. If K is 
slowly varying function and could be regarded to be a constant, the total F(R) behaves as 
F(R) ~ R m , which gives 

h = \, ^- (m - 1)(2 'T 1) . 

t m — 2 
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when hn > or 



t, — t m — 2 



when /i < 0. Here, h is a constant. 

Alternatively, the addition of an extra scalar field may bring the evolution to the phantom 
era. This is because such terms can become larger at the constant curvature and hence may 
induce such a crossing in realistic models. 

At present, it is not so clear if such non-local models could be variable or not due to 
technical problems, but this model has a possibility to explain the complicated cosmic ex- 
pansion, especially the coincidence problem. We now consider the models including such 
non-local terms to show the generality of our method. 

We remark that as shown in Ref. 26[, such a theory may successfully pass the solar 



251 ] . such non-local models may correctly 



system tests, and that as demonstrated in Ref. 
reproduce the whole expansion history of the universe expansion, at least in their equivalent 
scalar-tensor form. 



VI. STABILITY UNDER A QUANTUM CORRECTION 

In this section, we examine the stability for the obtained solutions of the crossing of 
the phantom divide under a quantum correction of massless conformally-invariant fields. It 
is convenient to do it by taking account of conformal anomaly induced effective pressure 
and energy-density. Note that we do not discuss the quantum regime of modified gravity 
itself because it is relevant at strong curvature (near to the Planck scale) where the form of 
modified gravity may be quite different from the one at the late universe. 

Quantum effects produce the conformal anomaly: 

T A = b (f + jjoR) + b'G + b"DR , (6.1) 

Here F : the square of 4d Weyl tensor, G : Gauss-Bonnet invariant, which are given as 

F = -R 2 — 2RijR l i + RijkiR 1 ^ 1 , G = R 2 — 4RijR 1 ^ + RijkiR^ kl . (6.2) 
o 

In the FRW background (12.61) . we find 

F = , G = 24 (HH 2 + . (6.3) 
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For N real scalar, N\/2 Dirac spinor, N\ vector fields, N2 (= or 1) gravitons and iVnD 
higher derivative conformal scalars, 

N + 6N 1/2 + 12JVi + 611iV 2 - 8A^ HD 
~ 120(4vr) 2 ' 

iV + lliVi/ 2 + 62JVi + 1411iV 2 - 28Nhd 

= - ; rr . 6.4 

360(4vr) 2 v ' 

b" can be arbitrary and we may choose, for example, b" = —26/3 or b" = 0. If we assume Ta 
can be given by the effective energy density pa and pressure pa from the conformal anomaly 

as 

T A = ~pA + 3pa , (6.5) 
and pa and pa satisfy the conservation law, 

Pa + 3H( P a+Pa) =0 , (6.6) 

we find 

PA = --A dta*HT A , p A = ~ [ dta A HT A + ^ . (6.7) 
Ob j «jd J o 

On the other hand, Eqs. (lC6jl and (1C7I) give the effective energy density pp and pressure 
Pf from f(R) = F(R) — R term: 

k 2 Pf = -- (F(R) -R) + 3(H 2 + H^j (F'(R) - 1) - 18 (aH 2 H + HH^j F"{R) , 
k 2 Pf = l - (F(R) -R)- (3H 2 + Hj (F'(R) - 1) 

+6 (8H 2 H + AH 2 + 6HH + Hj F"{R) + 36 UhH + &\ 2 F"\R) . (6.8) 

We now investigate the magnitude of pa, Pa, Pf, and pp when the phantom crossing 
occurs, when H = 0. We assume the magnitude of the Hubble rate H could be the order of 
the present Hubble constant H : 

H ~ H ~ 10~ 33 eV . (6.9) 

Expressions on (16. 7p tells that we may assume Pa ~ Pa ~ Ta- Then we find 

Pa~Pa~CH%. (6.10) 

Here C is a dimensionless constant coming from 6, b', b", and numerical constants and 
therefore C ~ 10 2 ~ 3 . On the other hand, expressions in (16.81) tell that we may assume 
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Pf ~ Pf ~ f(R)/ 1 ^ 2 - Since f(R) plays the role of the effective cosmological constant, we 
also assume f(R) ~ -ffg an d we fi n d 

p F ~ p F ~ — . (6.11) 

Since 1/k ~ 10 28 eV, we find 

\Pf\ > \pa\ , > \pa\ • (6.12) 

Therefore the quantum correction could be small when the phantom crossing occurs and the 
obtained solutions of the phantom crossing in this paper could be stable under the quantum 
correction. We should note that the quantum correction becomes important near the Big 
Rip singularity, where the curvature becomes very large. 



VII. CONCLUSION 



In the present paper, we have considered a crossing of the phantom divide in modified 
gravity. We have reconstructed an explicit model of modified gravity in which a crossing of 
the phantom divide can occur by using the reconstruction method proposed in Ref. 



. As a 



result, we have shown that the (finite-time) Big Rip singularity appears in the reconstructed 
model of modified gravity (i.e., in the Jordan frame), whereas that in the corresponding scalar 
field theory obtained through the conformal transformation (i.e., in the Einstein frame) the 
singularity becomes the infinite-time one. Furthermore, we have investigated the relations 
between scalar field theories with realizing a crossing of the phantom divide and the cor- 
responding modified gravitational ones by using the inverse conformal transformation of 
scalar field theories. It has been demonstrated that the scalar field theories describing the 
non-phantom phase (phantom one with the Big Rip singularity) can be represented as the 
theories of real (complex) F(R) gravity through the inverse (complex) conformal transfor- 
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which 



mation. Moreover, taking into account the fact that in the viable models 
are very realistic modified gravities that evade solar-system tests, our universe is asymp- 
totically de Sitter space, we have also proposed a model of modified gravity in which the 
transition from the de Sitter universe to the phantom phase can occur. We have found that 
to construct a viable model generating the transition from de Sitter space to the phantom 
phase, additional non-local term or almost equivalent scalar field is necessary. It would be 
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interesting to reconsider this problem in the presence of ideal fluid matter. In addition, we 
have examined the stability for the obtained solutions of the crossing of the phantom divide 
under a quantum correction coming from conformal anomaly. 

The study of the future evolution of the universe as discussed in this paper may be 
important to understand whether our universe evolves eternally or it will enter into the 
finite-time singularity. Additionally, this may shed extra light to some specific properties of 
different dark energy models and may help in selecting correct descriptions for the dark side 
of the universe. 

The originality of this work is to reconstruct an explicit model of modified gravity in 
which a crossing of the phantom divide can be realized. This is the point beyond the 
already existing literature. We have demonstrated that in principle the crossing of the 
phantom divide can occur at the present time or in the near future in the framework of 
modified gravity without introducing any extra scalar components with the wrong kinetic 
sign such as a phantom. This corresponds to the proof that the crossing of the phantom 
divide is possible also for modified gravity theories, which seem to be much less pathologic 
than the usual phantom scalar models of dark energy, similar to the scalar field theories in 
the framework of general relativity. The demonstration in this work can be regarded as a 
significant step to construct a more realistic model of modified gravity to correctly describe 
the expansion history of the universe. 
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APPENDIX A: NOTE ON THE RECONSTRUCTION METHOD 



In this appendix, we note the following point on the reconstruction method explained in 
Sec. II A. 

In the action (12.21) . if we redefine the auxiliary scalar field by <fi = by using a 

proper function $, and define P{p>) = P(&(ip)) and Q(p) = Q(§(p>)), the action 

S = j ^oV^ j_L [p(<p)R + Q(<p)\ +£ mattcr j (Al) 

is equivalent to the action (12.21) because this gives identical F(R) gravity. This can be 
explicitly confirmed as follows. First we should note that Eq. (I2.3P is modified as 

o = + = J dp (^)) R + fflil *t . (A2) 

dp dp \ rf$ d& f dp ' 

Then we can solve p with respect to R by 

up = cp(R) = . (A3) 

Here, p> is the inverse function of $. Hence the obtained F(R), corresponding to (12 .4p . is 

F(R) = P(p(R))R + Q(p(R)) = P ($ ($- 1 (0 (A))) i2 + Q ($ (^ >_1 (0 (i?))) 

= P((f>(R))R + Q(0(i2)) = F(i2) . (A4) 

Thus the obtained -F(-R) could be identical. Consequently, there are always ambiguities for 
the choice in <p like a gauge symmetry. In the FRW universe, we now assume that we have 
solved the F(i?)-gravity theory and obtained R as a function of time t as R = R(t). Then 
can be expressed as a function of t, = <f>(t), by a proper function <fi. If we redefine a scalar 
field hy p = _1 (0) by using the inverse function <p~ l of 0, we obtain p> = t. Hence we can 
always, at least locally, identify <p with time t, <p = t, which can be interpreted as a gauge 
condition corresponding to the reparameterization of <p = 4>(<p). 

There could be several cases that we cannot construct F(R). One possibility is that the 
differential equation (12. lip has no consistent solution. Another possibility could be the case 
that the algebraic equation (12.31) has no solution for obtained P and Q (for example, P 
and/or Q is a constant). 
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APPENDIX B: RECONSTRUCTION OF AN EXPLICIT MODEL 



In this appendix, we demonstrate that Eq. (12.131) can be a solution of Eq. (12. lip without 
matter. 

We start with Eq. (12. lip without matter: 

^p(0) d-g^dp^) 

= ^ — w^ +2 ^ 2 ~ m - (B1) 

By redefining P(<p) as Eq. (12.131) . Eq. (IB 1[) is rewritten to 

i ^n»*(*-™n , (B2) 



d(f) 2 

We now consider the model Eq. (12.141) . In this case, Eq. (1B2[) is reduced to 

1 d 2 p((p) 257(7 + 1) 



(B3) 



p(4>) d(f) 2 (f) 2 

which can be solved as Eq. (I2.15p . 

We mention the following point about the form of g(4>) in Eq. (I2.14p . From Eq. (12.191) . 
we see that as the universe evolves, the sign of H has to change in time so that a crossing 
of the phantom divide can occur. To realize such a behavior of H, there must exist (at 
least) two terms of in the brackets [] of the logarithmic function on the right-hand side 
of Eq. ( I2.14p . (Incidentally, the reason why we select the coefficient '10' on the right-hand 
side of Eq. (2.20) is to obtain the solution analytically.) As another form of the Hubble rate 
realizing a crossing of the phantom divide, there is the following model [ijj]: 

H = **{i + i^)' < B4) 

where ho is a positive constant. It has been shown that this cosmology can be constructed 



in terms of multiple scalar field theories 15]. This form also consists of two terms in t. In 
this model, H = h (2t - t s ) t s / [t 2 (t s - tf] . It follows from Eq. (12191) that when t < t s /2, 
H < and hence the universe is in non-phantom phase (w c g > —1), but that when t > t s /2, 
H > and thus the universe is in phantom phase (w c q < — 1). It is not hard to formulate 
the explicit model of modified gravity with the above type of the crossing of the phantom 
divide. 
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APPENDIX C: DERIVATION OF THE CORRESPONDING SCALAR FIELD 
THEORY 



In this appendix, we derive the expression of the action (13.41) . 

It follows from the action (13.1 p that the equation of motion of the other auxiliary field ( 
is given by 



(CI) 



where the prime denotes differentiation with respect to (. Substituting Eq. flCll) into 
Eq. (13.11) and eliminating £ from Eq. (13.11) . we find 

1 



S 



d x\J — g 



2k 



2 (F'(0R+F(0-F'(0() + £ mattcr 



(C2) 



We make the conformal transformation (13. 2p with Eq. (I3.3P of the action (1C2j) . Consequently, 



39, 



the action in the Einstein frame is given by 

S E = J d 4 xy^g~ 
where 

V{a) = e-°C{a) - e~ 2 °F (((a)) 



48] 



c 



p~ 2a r 

° '—matter 



(C3) 



no 



(C4) 



no (F'(o) 2 ' 

and g is the determinant of g^ u . In deriving Eqs. (1C3|) and (1C4I) . we have used Eq. (13.31) . In 
addition, ((a) in Eq. (IC4I) is obtained by solving Eq. (I3.3P with respect to £ as C = Civ 9 )- By 
defining y9 as tp = a/3/2ct/k, the action (1C3h is reduced to the form of the canonical scalar 
field theory M . 

From the action (12.11) . we find that the gravitational field equation is given by 

F'{R)R I1V - \g,„F{R) + g, v UF'{R) - V,V U F'{R) = k 2 T^) . (C5) 



When there is no matter, in the FRW background (12.61) the (fi, v) = (0, 0) component 
and the trace part of the (/i, v) = component of Eq. (lC5p . where i and j run from 1 to 
3, are given by 



3H 2 



- (F(R) -R) + 3[H 2 + H) (F'(R) - 1) - 18 ( AH 2 H + HH) F"(R) , (C6) 
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and 

- (2H + 3H 2 ^j = - (F(R) -R)- (3H 2 + Hj (F'(R) - 1) 

+ 6 (%H 2 H + AH 2 + 6HH + Hj F"(R) + 36 (aHH + H)* F"'(R) , (C7) 

respectively. Using Eqs. (1C6I) and (1C7I) . it follows 

2HF'(R) + 6 (-AH 2 H + AH 2 + 3M + jf ) ^"(.ff) + 36 (4M + Hj 2 F'"(R) = . (C8) 

APPENDIX D: CORRESPONDENCE BETWEEN THE SCALAR FIELD THE- 
ORIES AND F(R) GRAVITIES 

In this appendix, we explore the scalar field theories with realizing a crossing of the 
phantom divide in the Einstein frame and consider the behavior of the corresponding theories 
of F(R) gravity. 

1. Scalar field theories 

The action of scalar field theories in the Einstein frame is given by 

R 1 



S<s> = d xy'—g 



(Dl) 



where u ($) is a functions of the scalar field $ and W (<&) is the potential of $. 
In the FRW background ( 12.61) . the Einstein equations are given by 



3 2 • 
—H 2 = p$ , -H = p$ + p$ , (D2) 



where the energy density p$ of the scalar field $ and the pressure p$ of it are given by 

p$ = -oj ($) $ 2 + W ($) , p$ = \j ($) $ 2 - W ($) , (D3) 

respectively. Using equations in flD2j) and (|D3j) . we obtain 

c^($)$ 2 = -— (D4) 
W ($) = — f 3# 2 + if) . (D5) 
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It is t 

/ ($) as 



le interesting case that u ($) and W ($) are defined in terms of a single function 
49] 

2 dl ($) 



d$ ' 



3/ 2 ($) 



Thus we can find the solutions 



dl{$) 



(D6) 
(D7) 



$ = t, H = I(t). 

In what follows, we consider the case in which these solutions are satisfied. 
If we define a new scalar field x as 



(D8) 



X 



d$y/\u ($) 



(D9) 



the action (IDip can be rewritten to the form in Eq. (14. ip . where the sign in front of the 
kinetic term depends on that of u(<&). If the sign of to ($) is positive (negative), that of 
the kinetic term is — (+). In the non-phantom phase, the sign of the kinetic term is always 
— , and in the phantom one it is always +. In principle, it follows from Eq. (1D9|) that 
$ can be solved with respect to \ as $ = Hence, the potential W(x) is given by 

In the case of the model explained in Sec. II B, it follows from Eq. ( 12. 18ft that 7 (<&) is 
given by 

2 7 +l" 



T+(7+l)(f 



1 



2 7 +l 



(D10) 



From the solutions in (1D8p . we find that Eq. ( IDlOp gives 

27+1 



' T 



7 + (7 + l)(i) 



1 



2 7 +l 



(Dll) 



In deriving the expressions in (1D1QI) and (IDlip . we have used Eq. (12.171) . In this case, from 



26 



Eqs. (1D6|) . (lD7l) and (|D9]) we find 
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u;($) 



7-47(7 + 1) (g) -(7+l)(f 



2(27+1) 




(D12) 



(D13) 



x 



k 2 $ 2 (x) 

7 (30 7 - 1) + 



*(x) 



2 7+l 



64 7 ( 7 + 1) f *M) 



27+1 



+ (7 + 1) (30 7 + 31) 



2(27+1)' 



(D14) 



As shown in Sec. II B, when t < t c , H < (non-phantom phase) and it follows from 
Eq. ( 1D4I) that uj > 0. In the non-phantom phase, the sign of the kinetic term in the action 
(14. ip is — . On the other hand, when t > t c , H > (phantom phase) and from Eq. (1D4I) we 
see that uj < 0. In the phantom phase, the sign of the kinetic term in the action (14.11) is +. 
When t = t c , uj = and the transition from non-phantom phase to phantom one occurs 41]. 



2. Corresponding theories of F(R) gravity 

We investigate the behavior of the modified gravity (14.81) in the limit of t — > t s for the 
case in which the potential W (x{R)) is given by Eq. (1D14I) . If the scalar field x is real, it 
follows from Eq. (14.91) that the scalar curvature R is not always real. In order for R to be 
real, the following condition should be satisfied j^ : 

e'V^x (^k 2 W( X ) -iVGK^^j = e-V^x (^k 2 W ( X ) + iVG*^^ j . (D15) 

This condition is reduced to 

1 dW( X ) 



2\ -nt&nx -k X . (D16) 

W( X ) d X 

Except for the form of W(x) satisfying Eq. (1D16[) . R is complex if X is rea l The form of 
W(x) in Eq. (ID14I) cannot satisfy Eq. (ID16I) . This applies to the case in which the Hubble 
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rate is given by Eq. (1B4[) . If the scalar field x is pure imaginary and expressed as \ = irj, 
where t] is a real scalar field, and the potential W(x) contains only even power of x, Eq. (14. 9 p 
is rewritten to 



R = e^V 2 / 3 ^ 



dW(-rf\ 
d(-r] 2 ) 



(D17) 



which tells that R is real. In fact, however, that the action ( 14.11) with the + sign of the 
kinetic term is reduced to [47] 



S n 



d xy'—g 



A _ l^rfcr, - W(- V 2 ) 



(D18) 



This action corresponds to the non-phantom (canonical) theory. As a consequence, if the 
scalar field theory describes the phantom phase with the Big Rip singularity, the correspond- 
ing theory of F(R) gravity is usually complex. (Note that counterexample is known j^ . 
but even in this case, when t — > t s , the modified gravity becomes complex.) When t — > t s , 
from Eq. (1D14I) W(x) diverges and hence Fp(R) also diverges. 

In the Einstein frame, the scalar field couples with matter and therefore, the frame could 
be unphysical. The coupling changes the scale of time interval. We usually measure the time 
by using electromagnetism. In the limit of neglecting the local gravity, the measured time 
corresponds to the cosmological time in the Jordan frame but not in the Einstein frame. 



APPENDIX E: EXAMPLES OF VIABLE F(R) GRAVITY MODELS 

In this appendix, we show examples of viable F(R) gravity models. 
The modified part in F(R) in the action (12.11) can be separated as 



F(R) = R + f(R) 



(El) 



where f(R) is an arbitrary function of the scalar curvature R. The interesting viable model 



proposed in Ref. 



361 ] is given by 



f(R)=f m {R) 



M 2 c 2 (R/M 2 Y M 2 c 2 M 2 c 2 /c 3 



+ 



c 3 (R/M 2 ) p +1 c 3 c 3 (R/M 2 ) p + 1 



(E2) 



where c 2 and c 3 are dimensionless constants, p is a positive constant, and M denotes a mass 
scale. In this model, when R/M 2 — » oo, fns(R) ~ —M 2 c 2 /c 3 = const. 
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The generalizations of the above model which admits the unification of early-time infla- 
tion with late-time acceleration (with intermediate radiation/matter dominance) have been 
proposed in Ref. 37||. The example of such a theory is given by 



f(R) 



^tanh f 



-a 
- oij ( tanh 



bp (R — Rp 
2 

bi(R-Ri) 



\- tanh 
+ tanh 



bpRp 
2 

biRi 



(E3) 



where «o, «/, bp and bi are constant, and Rp and Rj are constant scalar curvatures. We note 
that such a theory may correctly describe the whole expansion history of the universe qual- 
itatively: inflation, radiation/matter dominance and dark energy, as it has been explained 



in Ref. 



39]. 



APPENDIX F: NUMERICAL CALCULATIONS WITH THE MATTER 



In this appendix, for the case Eq. (15. ip with t s gp — gi, we explain the numerical calcula- 
tions of Eq. (12. lip with the matter in detail. 

We define g = t s g and Y = 1 - X with X = t/t s . From Eqs. fF2TTUD and with 
t s g = 0i, H = dg(t)/ (dt) and R = 6 (h + 2# 2 ), we obtain 

(Fl) 
(F2) 

(F3) 



~g(t) = -~9o (Y + \ogY) 
a(t) 







90 










exp 


9o 









R = tiR = 6~g (2~g + ^ 



4gp , 2^0 + 1 



Y 2 



where we have taken a — (1 — l/a) go exp [g (1 — l/cc)] so that the present value of the scale 
factor should be unity. 

The solution of Eq. (1F3|) with respect to Y is given by 



Y(R) 



2~g 2 - R/6 



-20o 2 ± 



-20o 3 



0o (2g + l)R 
6 



(F4) 



In what follows, we use the lower sign in Eq. (JF4J). 

For simplicity, we consider the case in which there exists a matter with a constant EoS 
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parameter w = p/p. In this case, by using Eq. (1F4|) . Eqs. (12.111) and (12.121) are rewritten to 
[12g (2g Q Y + 2g + l)f 



l2 d 2 P(R) 



d 2 R 

~g Y 2 {12 [4~g Y + 3 (2~g + 1)] + Y 3 (1 + Y)} + 2~g Y*P(R) 

dR 



+ Y 6 t 2 s K 2 (l+w)p 



1 — —\ — 

a Y 



exp 



So 1 Y 

a 



-3(1+™) 



(F5) 



and 



t 2 s Q{R) 



1 + Y 



Y 



P(R) - 72~g 2 



1 + Y 
Y 3 



2^c 



2g Q + l\ dP(fl) 



dP 



+ 2t 2 s K 2 p 



'-a Y 



90 



exp 



~go 1 y 



-3(1+id) 



(F6) 



respectively. Here, p corresponds to the present energy density of the matter. In particular, 
we use the present value of the cold dark matter with w = for p, i.e., p = 0.233p c j^, 
where p c = 3H^/ (8nG) = 3.97 x 10~ 47 GeV 4 is the critical energy density. From Eq. (ED 



and R = t 2 R, we have 



F(R) 



(p(R)R + t 2 s Q(R)) ■ 



(F7) 



2k 2 2kH 2 s V 
To examine F(R), we numerically solve Eqs. (IF5I) - (IF7|) . 

In Fig. 3, we depict P (R) and Q(R)/v 2 as functions of R. The range of R is given by 
30 < R < 1000, corresponding to < X = t/t s < 1. Here, we have taken the initial 
conditions as P(R = 30) = 1.0 and dP{R = 30)/(dR) = 1.0. By using Eq. (|F7]), we show 
the behavior of F(R)/(2k 2 ) in Fig. 2. 



[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003); H. V. Peiris 
et al. [WMAP Collaboration], ibid. 148, 213 (2003); D. N. Spergel et al. [WMAP Collabora- 
tion], ibid. 170, 377 (2007). 

[2] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009). 

[3] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 
(1999); A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998); 
P. Astier et al. [The SNLS Collaboration], Astron. Astrophys. 447, 31 (2006); A. G. Riess et 
al, Astrophys. J. 659, 98 (2007). 

30 




FIG. 3: P(R) and Q{R)/v 2 as functions of R for t a = 2t and p = 0.233p c . 

[4] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003); V. Sahni, AIP Conf. Proc. 

782, 166 (2005) [J. Phys. Conf. Ser. 31, 115 (2006)]. 
[5] T. Padmanabhan, Phys. Rept. 380, 235 (2003). 

[6] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006). 
[7] S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007) 
[arXiv:hep-th/0601213] . 

[8] S. Nojiri and S. D. Odintsov, larXiv:0801. 48431 [astro-ph]; larXiv:0807.0685l [hep-th]; 

T. P. Sotiriou and V. Faraoni. larXiv:0805. 17261 [gr-qc]; F. S. N. Lobo. larXiv:0807.1640l [gr-qc]; 

S. Capozziello and M. Francaviglia, Gen. Rel. Grav. 40, 357 (2008). 
[9] U. Alam, V. Sahni and A. A. Starobinsky, JCAP 0406, 008 (2004); S. Nesseris and 

L. Perivolaropoulos, ibid. 0701, 018 (2007); P. U. Wu and H. W. Yu, Phys. Lett. B 643, 

315 (2006); H. K. Jassal, J. S. Bagla and T. Padmanabhan, |arXiv:astro-ph/0601389 
[10] R. R. Caldwell, Phys. Lett. B 545, 23 (2002); S. Nojiri and S. D. Odintsov, Phys. Lett. B 

562, 147 (2003) |arXiv:hep-th/0303117| . 
[11] L. Perivolaropoulos, JCAP 0510, 001 (2005); K. Nozari and S. D. Sadatian, Eur. Phys. J. C 

58, 499 (2008). 

[12] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2005) 



31 



[arXiv:hep-th/0504052] ; M. Sami, A. Toporensky, P. V. Tretjakov and S. Tsujikawa, Phys. 
Lett. B 619, 193 (2005); B. M. Leith and I. P. Neupane, JCAP 0705, 019 (2007); S. No- 
jiri, S. D. Odintsov and M. Sami, Phys. Rev. D 74, 046004 (2006) |arXiv:hep-th/0605039| ; 
T. Koivisto and D. F. Mota, Phys. Lett. B 644, 104 (2007); M. R. Setare and E. N. Saridakis, 
Phys. Lett. B 670, 1 (2008); A. K. Sanval. la7Xiv:0710.2450l [astro-ph], 
[13] A. Vikman, Phys. Rev. D 71, 023515 (2005). 

[14] M. z. Li, B. Feng and X. m. Zhang, JCAP 0512, 002 (2005) |arXiv:hep-ph /0503268| . 

[15] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71, 063004 (2005) 

[arXiv:hep-th/0501025] . 
[16] S. Nojiri and S. D. Odintsov, Phys. Rev. D 72, 023003 (2005) |arXiv:hep-th/0505215| . 
[17] L. P. Chimento and R. Lazkoz, Phys. Lett. B 639, 591 (2006); S. Sur and S. Das, JCAP 0901, 

007 (2009). 

[18] E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 043539 (2004) 
|arXiv:hep-th/0405034| ; S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 632, 597 
(2006) [arXiv:hep-th/0507182] ; E. Elizalde, S. Nojiri, S. D. Odintsov, D. Saez and V. Faraoni, 
Phys. Rev. D 77, 106005 (2008) |arXiv:0803.13TT1 [hep-th]]. 

[19] R. R. Caldwell and M. Doran, Phys. Rev. D 72, 043527 (2005). 

[20] H. Wei and R. G. Cai, Phys. Lett. B 634, 9 (2006). 

[21] B. Feng, X. L. Wang and X. M. Zhang, Phys. Lett. B 607, 35 (2005); Z. K. Guo, Y. S. Piao, 
X. M. Zhang and Y. Z. Zhang, ibid. 608, 177 (2005); X. F. Zhang, H. Li, Y. S. Piao and 
X. M. Zhang, Mod. Phys. Lett. A 21, 231 (2006); H. Mohseni Sadjadi and M. Alimohammadi, 
Phys. Rev. D 74, 043506 (2006); X. Zhang, ibid. 74, 103505 (2006). 

[22] W. Fang, W. Hu and A. Lewis, Phys. Rev. D 78, 087303 (2008). 

[23] M. Cataldo and L. P. Chimento. larXiv:0710.43"06l [astro-ph]. 

[24] B. Mclnnes, Nucl. Phys. B 718, 55 (2005); R. G. Cai, H. S. Zhang and A. Wang, Com- 
mun. Theor. Phys. 44, 948 (2005); R. G. Cai, Y. g. Gong and B. Wang, JCAP 0603, 006 
(2006); I. Y. Aref'eva, A. S. Koshelev and S. Y. Vernov, Phys. Rev. D 72, 064017 (2005); 
G. Kofinas, G. Panotopoulos and T. N. Tomaras, JHEP 0601, 107 (2006); I. Y. Aref'eva 
and A. S. Koshelev, ibid. 0702, 041 (2007); L. P. Chimento, R. Lazkoz, R. Maartens and 
I. Quiros, JCAP 0609, 004 (2006); P. S. Apostolopoulos and N. Tetradis, Phys. Rev. D 74, 
064021 (2006); Y. f. Cai, M. z. Li, J. X. Lu, Y. S. Piao, T. t. Qiu and X. m. Zhang, Phys. 



32 



Lett. B 651, 1 (2007); S. F. Wu, A. Chatrabhuti, G. H. Yang and P. M. Zhang, ibid. 659, 
45 (2008); K. Nozari and M. Pourghasemi, JCAP 0810, 044 (2008); K. Nozari, M. R. Setare, 
T. Azizi and N. Behrouz. larXiv:0810.142~7l [hep-th]: J. Sadeghi, M. R. Setare, A. Banijamali 
and F. Milani, Phys. Lett. B 662, 92 (2008). 

[25] S. Jhingan, S. Nojiri, S. D. Odintsov, M. Sami, I. Thongkool and S. Zerbini, Phys. Lett. B 
663, 424 (2008) |arXiv:0803.26l3l [hep-th]]. 

[26] S. Deser and R. P. Woodard, Phys. Rev. Lett. 99, 111301 (2007). 

[27] P. Singh, Class. Quant. Grav. 22, 4203 (2005); D. Samart and B. Gumjudpai, Phys. Rev. D 
76, 043514 (2007). 

[28] H. Stefancic, Phys. Rev. D 71, 084024 (2005); J. Grande, J. Sola and H. Stefancic, JCAP 

0608, 011 (2006). 
[29] M. Kunz and D. Sapone, Phys. Rev. D 74, 123503 (2006). 

[30] A. A. Andrianov, F. Cannata and A. Y. Kamenshchik, Phys. Rev. D 72, 043531 (2005). 
[31] M. C. B. Abdalla, S. Nojiri and S. D. Odintsov, Class. Quant. Grav. 22, L35 (2005) 

[arXiv:hep-th/0409177] . 
[32] L. Amendola and S. Tsujikawa, Phys. Lett. B 660, 125 (2008). 

[33] S. Nojiri and S. D. Odintsov, Phys. Rev. D 74, 086005 (2006) |arXiv:hep-th/0608008 |; J. 

Phys. Conf. Ser. 66, 012005 (2007) [arXiv:hep-th/0611071] ; J. Phys. A 40, 6725 (2007) 

[arXiv:hep-th/0610164] . 
[34] B. Mclnnes, JHEP 0208, 029 (2002). 

[35] R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 071301 (2003). 
[36] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007). 

[37] S. Nojiri and S. D. Odintsov, Phys. Lett. B 657, 238 (2007) }arXiv:0707.194Tl [hep-th]]: Phys. 
Rev. D 77, 026007 (2008) [arXiv:0710.1738l [hep-th]]; G. Cognola, E. Elizalde, S. Nojiri, 
S. D. Odintsov, L. Sebastiani and S. Zerbini, ibid. 77, 046009 (2008); |arXiv:0712.40T7l [hep- 
th]]. 

[38] S. A. Appleby and R. A. Battye, Phys. Lett. B 654, 7 (2007); L. Pogosian and A. Silvestri, 
Phys. Rev. D 77, 023503 (2008); S. Nojiri and S. D. Odintsov, Phys. Lett. B 652, 343 (2007) 
|arXiv:0706. 13781 [hep-th]]; S. Capozziello and S. Tsujikawa, Phys. Rev. D 77, 107501 (2008). 

[39] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003) [arXiv:hep-th/0307288l . 

[40] S. M. Carroll, V. Duvvuri, M. Trodden and M. S. Turner, Phys. Rev. D 70, 043528 (2004). 

33 



[41] S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 38, 1285 (2006) [arXiv:hep-th/0506212] . 

[42] S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 632, 597 (2006) 

|arXiv:hep-th/0507182] . 
[43] V. Sahni, Prog. Theor. Phys. Suppl. 172, 110 (2008). 

[44] S. Nojiri, S. D. Odintsov and M. Sami, Phys. Rev. D 74, 046004 (2006) 
|arXiv:hep-th/0605039] . 

[45] J. D. Barrow, Class. Quant. Grav. 21, L79 (2004); Phys. Lett. B 235, 40 (1990); S. Cotsakis 
and I. Klaoudatou, J. Geom. Phys. 55, 306 (2005); S. Nojiri and S. D. Odintsov, Phys. Rev. 
D 70, 103522 (2004) [arXiv:hep-th/0408170 ]; J. D. Barrow and C. G. Tsagas, Class. Quant. 
Grav. 22, 1563 (2005); M. P. Dabrowski, Phys. Rev. D 71, 103505 (2005); M. P. Dabrowski, 
Phys. Lett. B 625, 184 (2005); C. Cattoen and M. Visser, Class. Quant. Grav. 22, 4913 

(2005) ; L. Fernandez-Jambrina and R. Lazkoz, Phys. Rev. D 70, 121503 (2004); 74, 064030 

(2006) ; P. Tretyakov, A. Toporensky, Y. Shtanov and V. Sahni, Class. Quant. Grav. 23, 3259 
(2006); M. Bouhmadi-Lopez, P. F. Gonzalez-Diaz and P. Martin-Moruno, Phys. Lett. B 659, 
1 (2008); I. Brevik and O. Gorbunova, Eur. Phys. J. C 56, 425 (2008). 

[46] K. Bamba, S. Nojiri and S. D. Odintsov, JCAP 0810, 045 (2008) [arXiv: 0807. 25751 [hep-th]]; 

S. Nojiri and S. D. Odintsov, Phys. Rev. D 78, 046006 (2008) [arXiv:0804.35T9l [hep-th]]. 
[47] F. Briscese, E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Lett. B 646, 105 (2007) 

[arXiv:hep-th/0612220] . 

[48] Y. Fujii and K. Maeda, The Scalar-Tensor Theory of Gravitation (Cambridge University Press, 
Cambridge, United Kingdom, 2003); K. i. Maeda, Phys. Rev. D 39, 3159 (1989). 

[49] S. Capozziello, S. Nojiri and S. D. Odintsov, Phys. Lett. B 634, 93 (2006) 
[arXiv:hep-th/0512ll8] . 

[50] T. Chiba, Phys. Lett. B 575, 1 (2003). 



31 



